Charge ordering in the intermediate valence magnet YbPd 
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Abstract 

An x-ray diffraction study reveals the charge ordering structure in an intermediate valence mag- 
net YbPd with a CsCl-structure. The valence of the Yb ions forms an incommensurate structure, 
characterized by the wavevector (±0.07 ±0.07 1/2) below 130 K. At 105 K, the incommensurate 
structure turns into a commensurate structure, characterized by the wavevector (0 1/2). Based 
on the resonant x-ray diffraction spectra of the superlattice reflections, the valences of the Yb 
ions below 105 K are found to be 3+ and 2.6±. The origin of the long wavelength modulation is 
discussed with the aid of an Ising model having the second nearest neighbor interaction. 
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I. INTRODUCTION 



Electronic charges sometimes show intriguing behavior induced by their mutual interac- 
tions in condensed matter. Charge ordering is one of the most widely studied phenomena 
in transition metal oxides^, organic conductors 2 and rare earth compounds-^. Normally, 
charge ordering is seen in materials with low conductivity, and the phenomenon involves 
a sudden increase in the electrical resistivity. However, YbPd is a good metal in all tem- 
perature ranges 6 - with resistivity of ~ 10 2 fiQcm, and is proposed to have a charge order 
accompanied by a decrease in resistivity.- 

YbPd is an interesting compound in many ways. In the Yb-Pd phase diagram, this 
CsCl-structured intermetallic compound stands at the border of Yb 2+ and Yb 3+ , which are 
stabilized in the Yb rich side and the Pd rich side, respectively™. The valence of the Yb ions 
in YbPd measured by Ljjj-edge x-ray absorption spectroscopy was actually fractional 2.8+ 
across all temperature ranges 6 . In the fractional valence state, it shows magnetic ordering at 
Tn =1.9 K.- Magnetic ordering in the fractional valence state is not very common, and we 
were therefore interested in studying the origin of this rare phenomenon. Mossbauer mea- 
surements below 4.2 K show two kinds of Yb, one magnetic and the other non-magnetic, 
coexisting in equal proportions.- There are five phase transitions, at Ti=125 K, T2=105 K, 
Tn =1.9 K, Tmh =0.6 K and Tml =0.3 K.- While the low temperature transitions below 
2 K are known to be magnetic transitions, the other two transitions have not yet been 
explained. The high temperature transitions are accompanied by strong anomalies in the 
material's specific heat, thermal expansion 6 and elastic constants^ -. Although these ob- 
servations collectively suggest that either T x or T 2 is the charge ordering temperature, no 
superstructure has been reported to date. 

Charge ordering can be studied by ordinary x-ray diffraction because of the difference 
in the ionic radius of Yb 2+ and Yb 3+ , which produces a Pd displacement in the charge 
ordering phase. A more relevant probe for the spatial valence arrangement is resonant x- 
ray diffraction, which is a combination of an x-ray diffraction technique with spectroscopy. 
Using the L in absorption edge energy, it is possible to study the valence arrangements of 
4/ electron systems^. We have succeeded in clearly observing the superstructure induced 
by charge ordering in YbPd by means of these bulk sensitive techniques. 



2 



II. EXPERIMENT 



Single crystal samples were grown by the self-flux method.™ The as-grown samples have 
facets parallel to the {001}-planes and are cube-shaped. The typical size of these crystals 
is 1 mm 3 . Non- resonant x-ray diffraction measurements were performed with a four-circle 
x-ray diffractometer attached to a Mo K a x-ray generator. The incident x-ray beam was 
monochromatized with a bent graphite monochromator, and the scattered beam was then 
measured with a charge-coupled device (CCD) camera or a point detector. To find the 
superlattice reflections, x-ray photographs were taken with 18 keV intense synchrotron ra- 
diation at the BL-8B beamline of the Photon Factory, KEK, Japan. The charge order was 
examined by the resonant x-ray scattering technique at the Yb Lm absorption edge. The 
measurement was performed at the BL-4C beamline of the Photon Factory with a four-circle 
diffractometer. The sample temperature was controlled by using a closed cycle refrigerator 
for all measurements. 

III. RESULTS AND ANALYSIS 
A. Lattice parameters 

Figure [TJ^a) shows the 8-26 line profiles of (006) c Bragg reflection at 140 K (above T%, high 
temperature (HT) phase), 115 K (between T\ and T 2 , medium temperature (MT) phase) 
and 95 K (below T 2 , low temperature (LT) phase) measured with the Mo K a x-rays. Here, 
the suffix c denotes the index in the high-temperature cubic phase. In addition to the 
doublet caused by the K a \ and K a 2 x-rays, clear peak splitting with an intensity ratio of 1:2 
was found below T\. This result indicates that the highly symmetrical YbPd deforms into 
a tetragonal structure below Ti, and the single crystal turns into a multi-domain crystal. 
Panel (b) shows the lattice parameters a t and c t obtained from the peak positions of the 
(600) t and (006) t reflections, where the suffix t denotes the index in the tetragonal phase. 
Based on the peak profile, we found that T\ is between 130 K and 135 K. Above Ti, a t 
is equal to c t because the HT phase is cubic. Both of these transitions induce volume 
expansions with cooling. This behavior agrees well with the slope of the phase boundary in 
the pressure-temperature phase diagram shown in the inset of panel (b).- 

The isotropic atomic displacement parameters given by B = 87r(u 2 ), where u denotes the 
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FIG. 1: (a) 8 -26 line profile of (006) c Bragg reflection at 140 K, 115 K and 95 K measured with 
Mo K a x-rays. The error bars are shorter than the symbol size. The thick and thin gray curves 
for the 95 K profile show the results of the peak separation using a double Gaussian for the K a \ 2 
doublet. The dashed lines are intended as visual guides, (b) Temperature variation of the lattice 
parameters. Inset: Phase diagram reported in refU. 

atomic displacement from the equilibrium position, for Yb and Pd at 300 K were estimated 
from the (00Z) C intensities. The parameter values were 0.7 A 2 and 1.4 A 2 , respectively. 
Because the Lindemann melting criterion predicts a value of 0.8 A 2 for B at 300 K, the 
B for Pd is exceptionally large. The large B value for Pd may thus originate from the 
fluctuation of the Yb radius caused by the valence fluctuation. 
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FIG. 2: (Color online) Oscillation photographs around (136) c taken in (a) the HT phase (150 K) 
and (b) the LT phase (80 K). (c) Observed intensity at 95 K plotted against the calculated intensity 
for the best fitting structure, (d) Schematic view of the LT phase structure. 

B. Crystal structure and charge ordering in the LT phase 

The synchrotron diffraction experiment was performed to search for superlattice reflec- 
tions. Figure [2] shows the oscillation photographs around (136) c taken in (a) the HT phase, 
and (b) the LT phase. Along with the peak splitting caused by the cubic-tetragonal phase 
transition, the superlattice reflections characterized by the wavevector (00|) c were observed. 



The positions of the observed superlattice reflections were examined by using the four- 
circle diffractometer. The superlattice reflections in the LT phase were found to be at 
(00|)t away from the fundamental Bragg reflections and no reflections were seen at the 
(|00) t positions. The peak widths of the superlattice reflections were the same as those of 
the fundamental Bragg reflections, which means that the correlation of the superstructure 
reaches a long range. A strong (0 6.5) t reflection was observed, whereas (6 |) t reflection 
was very weak. This feature indicates that the atomic displacement 8 in the LT phase 
is parallel to the c-axis, because the superlattice reflection intensity caused by the small 
atomic displacement is proportional to \Q • S\ 2 , where Q denotes the scattering vector.— We 
measured the Bragg and superlattice intensities along the (00/) t line at 95 K to obtain the 
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structure in the LT phase. Figure [2^c) shows the integrated intensities of the (00Z) t peaks 
plotted against those fitted by the Bs and the <5s along the c-axis for both Pd and Yb. It 
was found that the Pd atoms are displaced by 0.1 lA along the c-direction to form a twofold 
structure in the LT phase, as shown in Fig. [2] (d). No atomic displacement was found for 
the Yb ions. This structure implies an alternate arrangement of high- and low-valence Yb 
ions along the c-direction. 

To clarify the valence arrangement, the energy spectra of the superlattice reflections were 
measured. Each element has a characteristic absorption edge, and the atomic form factor 
/ strongly depends on the x-ray energy E around the edge. Therefore, / is written as 
f(Q,E) = fo(Q) + f'(E) + if "(E), where f denotes the Thomson scattering term, and 
/' and f" denote the real part and the imaginary part of the anomalous dispersion term, 
respectively. The difference in the edge energies between the divalent and trivalent Yb ions 
is reported to be 7 eV.- This energy difference enables us to distinguish the valence values 
of these ions. 

The structure factor F for the superlattice reflections (00-|)t is written as 

i^oof = (-l) (i - 1)/2 2/ Pd sin(47rL5) + f H - f L , (1) 
f H = (n f + A)/3+ + {1 - (n f + A)}/ 2+ , (2) 
f L = (n f - A)/ 3+ + {1 - (n f - A)}/ 2+ , (3) 

where /pa, fz+ and f%+ denote the atomic form factors for the Pd, Yb 3+ and Yb 2+ ions, 
8 shows the Pd displacement in the unit of c t , rif is the average hole concentration of 
the Yb /-orbital with a value close to 0.8,- and A denotes the amplitude of the charge 
ordering. The energy dependence of the scattered intensity is \F\ 2 / n(E), where fi(E) denotes 
the absorption coefficient.— Using the theoretically calculated anomalous dispersion term 
for an isolated Yb atom^ (f c (E)) and the structure factor from eq.([T]), we analyzed our 
experimental results. The values of / 3+ and / 2+ are obtained from fo(Q) + f c (E ± 3.5eV), 
and n(E) was extracted from the fluorescence of the sample (see Figj3](a)). The theoretically 
calculated spectra were convoluted with the experimental resolution function, which was a 
Gaussian with a full-width at half maximum of 7 eV. The measured energy spectrum for 
the (00|)t superlattice reflection is shown in Fig. El^b). The upturn in the intensity with 
increasing energy toward the absorption edge is not expected without charge ordering, as 
indicated by the thin solid curve, and is therefore a signature of the charge order. Based 
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FIG. 3: (Color online) (a) Anomalous dispersion terms /' and /" for Yb 2+ and Yb 3+ based on 
Hartree-Fock calculations.— The linear absorption coefficient that was estimated from the fluores- 
cence spectrum is also shown, (b) Energy spectrum of the (00^) superlattice reflection measured 
at 80 K. The gray, thick solid, and thin solid curves show the calculated spectra for A=0.5, 0.2 
and 0, respectively. 

on a comparison between the experimentally observed spectra and the calculated spectra, 
we conclude that A ~ 0.2, i.e., that the high- and low- valence Yb ions are 3+ and 2.6+, 
respectively, for the LT phase of YbPd, as shown in Fig. [2] (d). 

C. Incommensurate structure in the MT phase 

In the MT phase, we discovered a splitting of the superlattice reflections. Figure @|a) 
shows the scattering intensity distribution of the (^2.5) t plane at 115 K. There are four 
peaks at (±x ±x 2.5) t . Panel (b) shows the line profiles along (£ £ 2.5) t (the dashed 
line in panel (a)) at 115 K (MT phase), and at 80 K (LT phase). The peak position, 
i.e., the incommensurability x, changes from 0.07 in the MT phase to with the phase 
transition at T 2 . This result indicates that the charge order in the MT phase is an alternating 
arrangement along the c-direction, which is similar to that of the LT phase, but it has an 




7 



in-plane modulation. Similar to the LT phase, the superlattice peaks are as narrow as those 
of the Bragg reflections, and therefore a long-range correlation was established. There are 
two possible modulation structures, depicted by the stripe and checkerboard structures that 
are presented in FigJSJ The stripe structure is characterized by valence modulation along 
the [110] t direction with a wavelength of ~ 7\/2a. This modulation structure provides two 
superlattice peaks at (+x +x 2.5) t and {—x —x 2.5) t in the region shown in Fig. SI (a). The 
other two peaks are produced by the twin variant with a modulation vector along [110] t . 
When we adopt the checkerboard structure, an alternating arrangement of 7a x 7a-sized 
tiles is expected to produce the four peaks shown in Fig. 0] (a). 

The temperature variations of the incommensurability and the intensity integrated over 
all the peaks around (0 2.5) t are presented in panel (c). The MT-phase superlattice 
intensity was found below 140 K, which is slightly above T\. The inset shows that the peak 
width of the superlattice intensity above T\ is broader than the instrumental resolution, 
meaning that the superstructure is short-ranged. The intensity above T\ is therefore given 
by the fluctuation. The intensity is proportional to the square of the Pd displacement, 
which must be a good indicator of the Yb valence. Below T2=105 K, the intensity is nearly 
constant, which means that the charge order amplitude in the LT phase is independent 
of temperature. The incommensurability varies gradually in the MT phase, and jumps to 
zero at T 2 . The absence of the higher harmonic peaks at (±m ±nx 2.5) t (where n is an 
integer) and the sudden decrease in the total intensity at T 2 with heating indicate that the 
modulation is sinusoidal. 



IV. DISCUSSION 



There are two possible interpretations of the intensity distribution in the MT phase: 
stripe and checkerboard structures. While we have no decisive data to distinguish the two 
structures, the stripe model is the more plausible of the two. In Fig. H] (a), there are some 
differences in the intensities of the four peaks. The intensity at (+x, +x, 2.5) t is equal to 
that at (—x, —x, 2.5) t , and the intensity at (+x, —x, 2.5) t is equal to that at (—x, +x, 2.5) t . 
However, the former pair is apparently stronger than the latter pair. This difference is 
accounted for by the domain ratio in the stripe model, while the checkerboard model cannot 
account for this intensity difference. 
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FIG. 4: (a) Scattering intensity distribution of the (£772. 5)t plane at 115 K. (b) Line profiles along 
the dashed line in panel (a) at 115 K and 80 K. (c) Temperature dependence of the incommensu- 
rability x and the integrated intensity. 

Next, we examine the origins of the incommensurate structure based on the stripe model. 
The valence states of the ions can be treated as Ising pseudospins. It has long been 
known that an Ising spin with a ferro-coupling for the nearest neighbor (Ji > 0) and an 
antiferro-coupling for the second neighbor (J 2 < 0), i.e., the axial next-nearest neighbor 
Ising (ANNNI) model, results in an intricate phase diagram that involves wide incommen- 
surate regions and commensurate regions separated by the first order phase transitions.— 1 ^ 
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FIG. 5: (Color online) Schematics of the possible modulation structures in the MT phase, (a) 
Stripe structure; (b) checkerboard structure. The dark and gray symbols represent the Yb ions 
with the different valence states. 



The ANNNI model has been applied to 4/ electron magnets including CeSbA£, CeBi^ and 
UPd 2 Si2^, to ferroelectrics such as NaNO^ and [N(CH 3 )4] 2 MCLi— ■ , and to charge order- 
ing systems such as NaV20^ and EuPtP 5 . Among these materials, EuPtP has significant 
similarities to YbPd, because both are metallic valence fluctuation compounds. 

Let us compare the known properties of the charge ordering in YbPd with those expected 
from the ANNNI model. First, the characteristic modulation vector varies from (0.07, 0.07, 
0) at Ti to (0.06, 0.06, 0) at T 2 . According to ref. [l5|, the wavelength is stabilized in a narrow 
temperature range with — J%j J\ ~ 0.27. When the temperature decreases, the calculation 
then predicts a first order phase transition to a ferroic arrangement within the c-plane, which 
was observed in our experiments at Ti- The ratio —J2/J1 is controlled by the application 
of pressure.— 1 ^ As shown in the inset of Fig. QJb), the LT phase vanishes around 2 GPa, 
whereas the MT phase remains above that pressure level. When we simply map the ANNNI 
phase diagram to the YbPd temperature-pressure phase diagram, then the charge order 
characterized by the wavevector (|||), which corresponds to 2v^2a in-plane periodicity, is 
expected for the ground state under pressures of more than 2 GPa. It is also expected 
that a large number of phase transitions will occur in the MT phase. We can test the 
applicability of the ANNNI model to YbPd by performing further diffraction measurements 
in high pressure environments. 

Another common origin of the long wavelength modulation is Fermi surface nesting. 
However, the nesting vector can vary only moderately. The cell volume change at T2 is 
almost the same as that within the MT phase, which means that the change in the Fermi 
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wavevector at T2 should be at a similar level to the change within the MT phase. As shown 
in the temperature dependence of the incommensurability x (Fig. 11(c)), the jump in x at T2 
is far greater than its shift in the MT phase. Therefore, the long wavelength modulation is 
unlikely to be caused by nesting of the Fermi surface. 

Finally, we discuss the effects of the charge ordering on the magnetic and conductive 
properties of the material. The application of pressure makes the LT phase unstable, and 
the ground state charge ordering structure is changed from an x = structure to an x 7^ 
structure (the MT structure) above 2 GPa. In accordance with this transition, Tn vanishes.— 
Therefore, the magnetic ordering at Tn requires a charge ordered structure in the LT phase. 
The charge ordering characterized by x = makes Yb 3+ and Yb 2,6+ sublattices, which can 
have different Kondo temperatures. The Yb 3+ sublattice can form a magnetic order because 
it is not in the mixed valence state. In the x 7^ structure in which the sinusoidal charge 
modulation develops, the Yb 3+ content is small, and therefore the magnetic ordering is 
suppressed. 

In general, charge ordering makes the carrier density decrease, and therefore the resistivity 
increases, as reported in Yb4As^ and in many charge ordered 3c? electron systems.- In 
contrast, the resistivity of YbPd decreases when the charge ordering occurs. This feature 
originates from the good metallic nature of YbPd. The carrier density is nearly unchanged 
by the charge ordering because the average valence remains constant, while the randomness 
of the valence arrangement, which scatters the carriers, is suppressed by the charge ordering. 

Summary 

We have performed a series of x-ray diffraction measurements on the valence fluctuat- 
ing compound YbPd. The material shows two-fold charge ordering characterized by the 
wavevector (00|)t below 105 K. Between 105 K and 130 K, the charge ordering structure 
is modulated, and the characteristic wavevector is (±x, ±2, |) t with x ~ 0.07. We propose 
that this long wavelength structure can be described by using the ANNNI model. 
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